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> 
O 

It- 

^—t , Aims. The effect of gravitational microlensing on the intensity of gravitational radiation as it propagates through an inhomogeneous 

\^ • medium is considered. Lensing by both stars and a power law spectrum of density perturbations is examined. 

Methods. The long wavelengths characteristic of gravitational radiation mandate a statistical, physical-optics approach to treat the 
effect of the lensing. 

Results. A model for the mass power spectrum of a starfield, including the effects of clustering and allowing for a distribution of stellar 
^ ■ masses, is constructed and used to determine both the amplitude of fluctuations in the gravitational wave strain and its associated 

^^f-^ ' temporal fluctuation spectrum. For a uniformly distributed starfield the intensity variance scales linearly with stellar density, cr, but is 

^■f^ I enhanced by a factor S crrf when clustering is important, where rp is the Fresnel scale. The effect of lensing by a power law mass 

1^-^ , spectrum, applicable to lensing by small scale fluctuations in gas and dark matter, is also considered. For power law mass density 

spectra with indices steeper than -2 the wave amplitude exhibits rms fluctuations 1.3(AS^)''''(/)eft71 Gpc)''-%, where (Air) is the 
variance in the mass surface density measured in pc"* and D^g is the effective distance to the lensing medium. For shallower spectra 
' the amplitude of the fluctuations depends additionally on the inner length scale and power law index of the density fluctuations. The 

C~>) , intensity fluctuations are dominated by temporal fluctuations on long timescales. For lensing material moving at a speed v across the 

— , ■ line of sight the fluctuation timescale exceeds u"'(£)cff/l)''-. Lensing by small scale structure induces at most » 15 % rms variations if 

i-S^ ' the line of sight to a gravitational wave source intersects a region with densities ~ lOOMopc"^, which are typically encountered in the 

Oh' vicinity of galaxy clusters. 

O ' Key words. Gravitational Lensing - Gravitational Waves - Scattering - Galaxies: structure - dark matter 

^' 

' 1. Introduction 

> 

• I— I , The amplification of electromagnetic radiation due to gravitational microlensing (Byalko 1969; Chang & Refsdal 1979, 1984; 
Paczyriski 1986) is now a well-established phenomenon that has proven to be a successful tool in probing planetary- to solar-mass 
objects in the vicinity of the Milky Way (e.g. Afonso et al. 2000). Gravitational radiation is also subject to lensing by small- 
. . . scale structure, but such lensing differs from its electromagnetic analogue in several qualitative respects. In particular, under most 
circumstances it demands a treatment based on physical optics rather than geometric optics, and the number of objects contributing 
instantaneously to the lensed signal can greatly exceed that expected in the electromagnetic case. 

The lensing of gravitational radiation occurs in a regime in which wave effects cannot be ignored (see Macquart 2004, hereafter 
M04, and references therein). Gravitational wave detectors such as LIGO operate at frequencies 10 - lO'^Hz, while the planned 
Laser Interferometer Space Antenna (LISA) will operate from 10 ' - 10 "^Hz. Thus the size of most objects likely to lens such 
radiation is large compared to the wavelength of the radiation. A sufficient (but not necessary) condition for diffractive effects to be 
important is that the wavelength exceeds the sizes of the lensing objects, which occurs for objects of mass $ 10^ v"' M© (Takahashi 
& Nakamura 2003). 

The radius of the region that contributes to wave fluctuations is characterised by the Fresnel scale, rp ~ '\jD/k, where D 
is an effective distance to the lensing material and k is the wavenumber (M04). This quantity exceeds 1 pc for typical radiation 
wavelengths and lensing geometries (i.e. with D characteristic of lensing at cosmological distances), so that many tens to thousands 
of objects (e.g. stars), depending on the wavelength, may contribute to the lensing. It is often therefore impractical to specify the 
exact mass distribution at all locations and all times, and the lensing of gravitational radiation lends itself to a treatment in which 
the mass density fluctuations of the lensing objects are described statistically. 

Most operational and proposed detectors are directly sensitive to the wave amplitude of the gravitational radiation itself, and 
many gravitational lensing effects relevant to these detectors are described in terms of second-order moments of the wavefield or 
mutual coherence functions. These were investigated in M04 (see also Takahashi & Nakamura 2003), where a statistical treatment 
was introduced to determine the effect of lensing on the properties of the gravitational radiation. 
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In the present paper this formalism is extended to compute the intensity fluctuations induced by gravitational lensing by small- 
scale structures. More generally, we consider the efl'ect of lensing on the fourth-order moment of the wavefield, of which a special 
case describes intensity fluctuations. There are several reasons to consider this quantity. 

The first is that the wave power, measured by the intensity, is a fundamental quantity that may be affected by gravitational 
lensing. Holz & Hughes (2005) suggest that binary black-hole mergers are suitable as "gravitational-wave cosmological sirens" 
if their electromagnetic counterparts can be identified because the frequency of the gravitational radiation and its time derivative 
in such systems uniquely determines the intrinsic luminosity. Uncertainties in the identification of the electromagnetic counterpart 
are estimated to introduce a 1-10% uncertainty in the distance relation (Holz & Hughes 2005). However, lensing-induced intensity 
variabiUty further hinders the suitability of such sources as standard sirens. Lensing potentially causes misestimation of the source 
intensity and, in long lived sources, may hinder identification of those variations which are intrinsic to the source itself. The effect 
of lensing by large-scale structures can be accounted for by examining the line of sight for lensing objects, but it is estimated that 
at least half the uncertainty is from lensing by structures less than ~ 0.6 Mpc in size (Kocsis et al. 2006). This subject has already 
received some attention through examination of fluctuations in the wave amplitude and phase (Takahashi 2006). 

The second is that the calculation of the ensemble-average mutual coherence presented in M04 is not sufficient to fully describe 
the effects imposed by gravitational lensing on the wave properties. Practical observations of gravitational wave sources extend only 
over a finite duration, so that the measured mutual coherence may deviate significantly from its average value. Thus comparison 
of the ensemble-average mutual coherence against its observed value may be invalid. The fourth order moment is an estimator 
of the amount by which the second moment instantaneously varies from its ensemble average value because, in its most general 
form, it may be regarded as the cross-correlation of two instantaneous mutual coherence measurements. We note that the fourth 
order moment is computed in other branches of scattering physics for this precisely this reason. This is particularly the case in 
scintillation physics where, even though a large ensemble of phase fluctuations may contribute instantaneously to the scattering, the 
instantaneous mutual coherence can exhibit large modulations around its ensemble average value (e.g. Goodman & Narayan 1989). 

A third is that lensing-induced intensity variations can potentially be used to recover the properties of the intervening lensing 
material. It was shown in M04 that the ensemble-average mutual coherence is related to the auto-covariance of lensing-induced 
phase delays, which is in turn related to the underlying mass power spectrum by a Fourier transform. It is pertinent to investigate 
whether fourth-order moments of the wave amplitude provide supplementary information on the underlying mass distribution. 

The physics of gravitational wave microlensing strongly resembles that encountered in the scattering of radio emission from 
pulsars and compact quasars as it propagates through density inhomogeneities in the intersteUar plasma (Rickett 1977). In light 
of this similarity, the terms microlensing and scattering are used interchangeably throughout the text. The main distinguishing 
characteristic of the present situation is that the power spectrum of lensing-induced phase fluctuations is typically steeper than that 
experienced in other applications, due to the manner in which the gravitational phase delays are generated from the underlying mass 
fluctuations. This qualitatively alters the nature of the intensity fluctuations. 

The large wavelengths typical of gravitational radiation, in the range lO^m ^ A 10^' m, introduce a number of physical 
efl^ects that are usually unimportant in the microlensing of electromagnetic radiation. In particular, the coherence area over which 
objects (e.g. stars) may contribute to the gravitational wave amphtude received at Earth can exceed many square parsecs so that, 
instead of a single object dominating the lensing at any given instant, an enormous number of objects can contribute to the lensing 
instantaneously. Thus we need to consider the simultaneous cumulative effect of a large number of objects on the wave amplitude. 

An important qualitative manner in which the lensing characteristics may change is through the effect of clustering. Since the 
efl'ective area that contributes to the wavefield encompasses many stars, stellar clustering becomes important in a way that is not 
necessarily the case in electromagnetic lensing. If the coherence area encompasses a cluster of stars we might expect the intensity 
deviations will more closely resemble lensing by a macrostar of mass NM whereas the incoherent contribution of individual stars 
randomly distributed in space might be only N^'^M if clustering were unimportant. 

The properties of the sources themselves also influence the nature of the gravitational lensing. In contrast to most lensed sources 
of electromagnetic radiation, sources of gravitational radiation are extremely compact and radiate "coherently", in the sense that 
they may be regarded as single emitting particles. Their compact nature is of particular importance, since it renders their radiation 
far more susceptible to certain lensing effects that are not relevant to large lensed sources. For large sources ampUfication effects 
associated with structure in the lensing pattern on angular scales less than the source angular diameter are heavily suppressed (Little 
& Hewish 1966; Narayan 1992); in effect the source "washes-out" lensing structure smaller than its own angular size. 

In this paper we derive the relation between the properties of the lensed gravitational radiation and the underlying mass distri- 
bution and use it to consider lensing by stars and a power law spectrum of mass density fluctuations. The layout of this paper is as 
follows. In section 2 we review the physics of gravitational wave propagation and employ it in section 3 to develop a statistical the- 
ory for the intensity fluctuations induced by lensing on a gravitational wave. This section also considers more general fourth-order 
moments of the wavefield, used to describe how the mutual coherence deviates from its mean quantity. The theory is then applied 
in Sect. 4 to lensing by a power law mass power spectrum and used to develop an intuitive understanding of the effect of lensing on 
gravitational wave properties. In Sect. 5 a model for the mass surface density power spectrum of a starfield is constructed, taking 
into account the clustering properties of the stars. This model is applied in Sect. 6 to derive the characteristics of the intensity fluc- 
tuations due to lensing by stars over a distribution of masses, and whose clustering properties may vary as a function of mass. The 
implications for the properties of gravitational radiation are presented in Sect. 7. The results are summarized in Sect. 8. Although the 
results are derived here in the context of gravitational wave lensing they are equaUy applicable to any situation in which the lensing 
field is best described using a statistical approach since, to lowest order in the gravitational wave amplitude, the propagation of both 
gravitational and electromagnetic radiation are described by identical wave equations. The physical optics treatment necessitated by 
the long wavelengths typical of gravitational radiation is, of course, also valid in the treatment of electromagnetic radiation. 
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2. Wave Propagation 

The amplitude of gravitational radiation propagates according to the following scalar wave equation (see e.g. Takahashi & Nakamura 
2003, M04), 

(V^ + a?)^ = Aa?U{r)^, (1) 

where u) - 2nv, f/(r) <k 1 is the gravitational potential, and 0(v, r) is the temporal Fourier transform of the scalar wave amplitude. 
We consider the thin lens geometry illustrated in Fig. 1, for which the solution of eq. ([TJ is (e.g. Schneider 1987) 

^(v, X) = — — y I cfn exp 



Inrl 

F 



I 



(x-X)2 + /(A(x) 

2rr 



F 



(2) 



where it is assumed that the radiation incident on the lensing plane emanates from a point source of unit amplitude. The scaled 
co-ordinate, X, is related to the position on the observer's plane, x, by X = Dis /Dsx, where Dls is the angular diameter distance 
from the lens to the source, and Dg is the angular diameter distance from the observer to the source. The angular diameter distance 
from the observer to the lens is Dl and zl is the associated redshift. Throughout the text all displacements on the observer's plane 
are written in terms of the scaled co-ordinate X. The Fresnel scale, defined as 

- ^ " IT.) 

is interpreted as the scale length on the lensing plane over which geometric phase delays (i.e. phase delays in the absence of mass 
perturbations) become important. The phase delay driven by fluctuations in the gravitational potential is 



<A(x) = K 



J^V2(x')ln(^), /r = -8.i^§, (4) 

where 2;(r) is the mass surface density projected onto the lensing plane. The scale factor xq is set to unity hereafter since only phase 
differences across the lensing plane cause effects of interest to us here, and the total phase delay is unimportant. 

The thin-lens approximation commonly employed in applications of electromagnetic lensing is still valid in the context of 
gravitational wave lensing despite that fact that the depths of the lensing objects may be much less than one wavelength. This is 
because the extremely small amplitude of the perturbation associated with a gravitational wave renders nonlinear effects associated 
with sharp gradients in the gravitational potential negligible (see, e.g., Thorne 1983). Although the lensing material may, in general, 
be distributed along the entire line of sight from the source to the observer, the thin-lensing approximation adopted throughout this 
paper is also sufficient for the lensing geometries considered here. This approach avoids the complexity introduced by the extra 
mathematical machinery needed to treat phase fluctuations caused by mass inhomogeneities distributed along the entire ray path. 
Moreover, extended medium treatments (Tatarski 1967; Codona & Frehlich 1987) show that in all cases of interest here it is possible 
to represent the scattering physics of an extended medium in terms of an equivalent thin screen (Tatarski & Zavorotnyi 1980). 



2.1. A statistical description of tine gravitational pinase delay 

An investigation of wave amplitude variations due to lensing by small scale structure requires a statistical description of the mass 
fluctuations that drive the phase perturbations. Consider matter distributed with volume density p(R) with corresponding power 
spectrum ^p{qx, qv, qz', z), where the direction of propagation is parallel to the z-axis, and it is assumed that the outer scale of the 
power spectrum is small compared to the total path length. The Markov approximation is assumed so that changes in the shape or 
amplitude of the power spectrum are parameterized as a function of distance along the propagation axis. The power spectrum of the 
mass perturbations projected onto an equivalent thin lensing plane is. 



(Ds(q) = J dz%(q,q,^Q;z), (5) 

where q is the reciprocal co-ordinate to position on the two-dimensional lensing plane. We make the simplifying assumption that 
the properties of the matter distribution on small scales are decoupled from the large-scale (i.e. galactic-scale) matter distribution, so 
that the mass surface density fluctuations are assumed to be wide-sense stationary. This approximation is valid for the small scales 
over which we wish to characterize the mass power spectrum. As such we explicitly ignore large-scale variations in the phase delay 
due to macrolensing by large objects, such as galaxies or clusters. Their contribution is not relevant to the physics considered here 
and their effect is easily introduced separately if necessary. 

The phase statistics are described by the phase autocovariance function, which quantizes the correlation of the phase between 
two points separated by a displacement r on the lensing plane. This correlation is related to the underlying power spectrum of mass 
surface density fluctuations by. 



Q(r) = mr' + r) - 0-)(^(r') - 0-)) ^ j d^q-^e-''^' (D£(q). 



(6) 



An equivalent but often more convenient quantity is the phase structure function, which measures the mean square phase differ- 
ence between two points. 



D^(r) = <[0(r' + r) - 0(r')]') = 2K^ J d\q-^ [l - e '"^'] ^^(q). 



(7) 
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There is a simple relationship between the phase structure function and the mean correlation between the wave amplitude measured 
between two receivers separated by a displacement r when the phase fluctuations obey Gaussian statistics (M04): 



<y(r)) = exp 



(8) 



where the incident radiation is again assumed to emanate from a point source of unit intensity. Eq. ([8]) shows that the mean intensity 
of the gravitational radiation, obtained by setting r = 0, is equal to the intrinsic source intensity. It also demonstrates that the average 
correlation between the wave amplitudes measured a vector r2 5^ apart is reduced by the phase and amplitude fluctuations induced 
by lensing. The value of this decorrelation contains information about the medium through which the radiation has propagated, and 
can potentially even be used to extract the power spectrum of mass fluctuations responsible for the lensing. However, eq. dHJ above 
represents an average over the ensemble of all possible lensing fluctuations consistent with a given mass power spectrum, O^- In 
a realistic situation in which the intensity or mutual coherence is averaged over a finite time interval, these quantities will deviate 
from their ensemble-average values. The following section considers the magnitude of these deviations and the timescales on which 
they occur. 



3. Intensity variance and other fourth moments of the wavef ield 

In this section we derive expressions for both the intensity variance and the power spectrum of temporal intensity fluctuations. The 
latter quantifies the timescale on which fluctuations are expected. We also consider deviations in the mutual coherence function, 
y(r), from its mean (ensemble-average) value by calculating both the variance and temporal power spectrum of its fluctuations. All 
quantities of interest here are computed by considering the following fourth-order moment of the wave amplitude, 

r4(ri,r2) = <M(r)M*(r + r2)M(r + ri)M*(r-i-ri -i- r2)). (9) 

This measures the correlation between the wave amplitudes measured by four hypothetical detectors whose locations are depicted 
in Figure 1. Various specializations of this quantity yield measures of particular interest to us. The choice ri = r2 = gives (/^), 
while allowing ri yields the correlation in the intensity, (/(ri + r')/(r')), between two hypothetical detectors separated by a 
displacement ri . Temporal intensity variations are calculated by making the replacement ri = vf. The choice ri = 0, r2 9^ yields 
the variance in the mutual coherence V(r2)- In the most general case, ri 0, r2 0, the quantity r4 measures the correlation in the 
mutual coherence between two pairs of receivers separated by a distance ri = vf, with each pair measuring a visibility V{r2)- 

To compute r4 we first consider the most general fourth moment of the wave field, which measures the correlation in the fields 
measured at the four points Xi, X2, X3 and X4 on the observer's plane. 



r(Xi,X2,X3,X4) = (0(Xi)0*(X2)0(X3)0*(X4)) = — ■ f c/Vjt/^Xjt/^Xjt/^x; eXp i [(x'l - Xi)^ - (x^ - X2)H (x', - X3)2 - (x^ - X4)^l 

+1x1 j d^I.(e [in II - x'll - In II - x^l + In II - x;,| - In || - x;|] 



Making the co-ordinate transformation, 

R' = ^[x'l -i- X2 -H X3 -H X4] 

r'l = i[x'i +x'2-x'^-x'^] 

p 



(11) 



with a corresponding change of variables for the unprimed coordinates, the fourth moment then becomes, 
r(xi,X2,X3,X4) = — ^ r dh\dh'2d^R'd^p'exp\-^[ir\-ri)-{r2-r'2) + ip-p')-(R-R')]] 

(exp iK j d|m)(ln||-x;|-ln||-x^| + ln||-x;,|-ln||-x;|) ^. 

(12) 

Cast in this form, we apply the simplifying assumption that the gravitational phase, i//, is wide-sense stationary, so that the statistics 
of the phase diff'erences are independent of the centre co-ordinate, R'. The integration over R' yields a factor e'^^P P^6(p - p'). 
Figure 1 shows that the four receivers are arranged in a parallelogram, implying p = 0, and integration over p' yields a factor 
{2nrh^ leaving. 



r4(ri,r2) = -r^-j-r { d^'^d^'^Qxp 



4-(r'i - ri) ■ (r2 - r'^) 



x(exp[i/:Z(r;,r^)]>, (13) 
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source 



lens plane velocity V^ff 



observer's plane 



lens plane 



Fig. 1. The fourth-order moment of the wave field, r4(r i , r2) measures the correlation between the two visibilities, each of length r2, 
measured a distance ri apart. Of the four hypothetical dectors depicted, detectors 1 and 1' measure the first visibility, while dectors 
2 and 2' measure the second. The quantity r4(0, ri) measures the variance in the visibility V(r2)- When the baseline r2 is set to zero, 
so that dector 1 is co-located with dector 1', and dector 2 is co-located with dector 2', the visibility V{r2) collapses to just the wave 
intensity. The quantity r4(ri,0) thus measures the autocorrelation of the intensity (/(r + ri)/(r)). The quantity r4(0, 0), when all 
four detectors are at the same location, measures {/^). 

where the effect of the gravitational phase delays is embodied in the quantity. 



Z(r;,r^) = Jd^lnl^l + ^i^j - + ^l^j + - ^i^j - - 



(14) 



In the following two subsections averages over the phase fluctuations are computed to derive expressions for the variance and power 
spectrum of both intensity and mutual coherence deviations. 



3.1. Small phase perturbations 

When the lensing-induced phase fluctuations are small the fourth moment is well-approximated by expanding eq. ( fT3] ) in powers of 
^Z(r;,r^), 



r4(ri,r2) = f dh\dh'2exp{-^ [{r\ - n) ■ (r2 - r^)]| [l + iK{Z(r[,r'^)} - ^{Zir'^y^fy+oU') 

{2nri,Y J [ri, j [ 2 



(15) 



This approximation is valid when the phase delay term KZ(r'^ , r'^) is less than unity over the region of integration. The eff'ective 
integration area is determined by the region over which the kernel exp[(ri - r[) ■ (r2 - r'2)/rp] varies slowly. This condition is 
satisfied when the rms phase delay across a distance of order the Fresnel scale on the lensing plane is less than unity. We shall see 
later that for most cases of interest the weak scattering approximation applies under far less restrictive conditions. 

No assumption is made about the statistics of the mass fluctuations except that they are wide-sense stationary and thus that the 
average mass density is independent of position, so one has (2(x)) = (2(x + x')) and (Z) = 0. It is easily shown that the expectation 
of Z^, written in terms of mass surface density autocorrelation functions, is 



<Z'(r,,r2)) 



dQ.dj3 In lal In \j3\ 



4Cda -13)- 2C^(a - /3 - r^) - 2C^(a - /3 + r^) - 2C^(a -f3-r\)- 2C^(a -f3 + r\) 



+CY(a. - /3 + r'l + r'2) + CY.(a - (B + r\ - r'2) + C-z(a - f3-r[+ Fj) + C^{a. - f3-r[- Tj) 



(16) 



The relationship between the mass density autocorrelation function and power spectrum in equation ^ allows (Z^) to be recast in 
terms of the power spectrum of mass density fluctuations. Subsequent integration over r'^ and yields our main expression for the 
fourth moment in the regime of weak scattering. 



r4(ri , rz) = J ?"'^%(q)(2 - 2 cos q ■ r2 - 4e"' '"' sin- 



q ■ (r2 - qrl) 



(17) 
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Equivalent results have been derived in other fields in which similar physics applies; it has been investigated extensively by Cronyn 
(1972) in the field of interplanetary scintillation. The main difference is that in the present case the phase delays are driven by 
distortions imposed by matter on the geometry of spacetime, whereas for interplanetary scintillation phase delays are related to 
density fluctuations by the plasma dispersion relation. 

Equation (fTTT i is the starting point for computing all quantities of interest in the weak scattering limit. In this section expressions 
for all these quantities are listed, but their explicit evaluation for given mass distributions is deferred to the next section. The most 
obvious quantity is the total variance of the intensity fluctuations. Since r4(0, 0) = (/^) the intensity variance is 



{An = r4(0,0)- 1 =4/:^ d\q-*(^^iq)sm 



2^2 \ 



(18) 



The timescale on which the intensity fluctuations occur is described by the power spectrum of the intensity fluctuations. It was 
mentioned above that the fourth moment gives the autocorrelation of the intensity fluctuations between points separated by a vector 
ri. The Fourier transform of this moment is used to obtain the power spectrum of intensity fluctuations across the observer's plane: 



W/(q) 



d ri exp[-/q ■ ri]r4(ri,0). 



(19) 



The integral of W7(q) over all q of course yields the variance in the intensity. In the weak scattering limit the spatial power spectrum 
of intensity fluctuations is derived from the fourth order moment by noting that eq. ( [TtT i is already written in terms of a Fourier trans- 
form. Setting r2 = and Fourier transforming r4(ri, 0) with respect to ri we see that the power spectrum of intensity fluctuations 



IS, 



W,(q) 



: AK^q-'^O^iq) sin^ 



2 2 



(20) 



Wave effects become particularly important for wavenumbers q greater than rp', where the power spectrum oscillates with a fre- 
quency determined by the Fresnel scale. On the other hand, the power spectrum is an unweighted measure of the power spectrum 



< r- 



In this 



of mass surface density fluctuations for wavenumbers q less than . In particular, one has W/(q) ~ K-'r^(^^{q) for q 
regime the power spectrum is also independent of observing frequency, with the frequency dependence of the Fresnel scale being 
offset by the dependence of K. At higher observing frequencies wave effects become less important, the Fresnel scale becomes 
smaller, and the regime of applicability of the approximation q < r^^ means that wave effects are important only for progressively 
smaller length scales. 

In practice it is difficult to measure the power spectrum of intensity fluctuations across a two-dimensional observing plane. A 
more practical observable quantity is the power spectrum of intensity variations experienced by a single detector. In general the 
lensing material (e.g. stars, gas and dark matter) will be in motion relative to the line of sight to the source. If the motions of the 
individual lensing constituents are assumed to be negligible relative to the bulk velocity of the lensing plane across the line of sight, 
the so-called frozen-screen approximation, the effective lensing velocity is. 



VEa.h|l-f^)+Vs.c(^ 



(21) 



where VEarth, Vscreen and Vsic are, respectively, the velocities of the Earth, lensing plane and the source transverse to the line of sight. 
The power spectrum of temporal intensity fluctuations is obtained from the fourth-order moment by equating spatial and temporal 
displacements according to the rule r — > vf to yielcQ, 



Wi(oj) 



dte"^' T^{\t, 0). 



(22) 



It is often the case that the power spectrum of mass fluctuations is an isotropic function of q, so that the fourth moment is also 
an isotropic function of ri. Without loss of generality the velocity can then be oriented along the jc-axis v - {v, 0), and the power 
spectrum of temporal intensity fluctuations is thus. 



Wiioj) 



dq 



-2 



(23) 



Both the variance and power spectrum of the mutual coherent fluctuations are derived in a similar manner to the intensity 
fluctuations. In subtracting out the mean mutual coherence to determine the variance, it is useful to recognise that the first three 
terms in eq. ( fTTT i just represent the square mean mutual coherence V(r2) in the small perturbation limit: 



{V(r2)f = exp[-D^(r2)/2]2 ^1-2K- 



d\q-*(:>^(q)[l - cos(r2 ■ q)] + 0(K*). 



(24) 



' Since coordinates on the observing plane are scaled (see the text below eq.|2}, <^ is also scaled. The replacement ai — > (Du /Ds )iD converts to 
observational quantities. This replacement is insubstantial in most cases under consideration since the modifying factor is of order unity. 
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The exponential function is expanded in the same way that it was linearised when making the weak scattering approximation for 
the fourth order moment. Making the separation - (V + 6V)^ - + V^, the fluctuation in the mutual coherence on a baseline 
r2 is AV\0, r2) = r(0, rz) - {Vf, 



j d\q-' 



<I)5;(q) sin^ 



q ■ (r2 - q'-p) 



(25) 



The power spectrum of mutual coherence fluctuations measured on a baseline ¥2 is 



dq 



— + r 



-2 



sm 



^r2, + qr2y-{^+q^)rl 



(26) 



where the mass fluctuation spectrum is again assumed to be an isotropic function of q. 



3.2. Large phase perturbations 



The regime of strong scattering applies when KZ(r'^ ^^2) ~ 1 o^^r '^he eff'ective region of integration in r'j and rl, . In order to compute 
averages over the phase in this regime, slightly more restrictive assumptions must be made about the statistical properties of the 
mass density fluctuations than is necessary in the regime of weak scattering. 

The gravitational phase delay is assumed here to be normally distributed. This is the case if the underlying mass surface density 
fluctuations are themselves normally distributed. Averages over the phase fluctuations are performed using the fact that, for a 
Gaussian random variable, 6x, with zero mean, the average of {e\p[i6x]} is e\p[-{6x^}/2]. This implies that the average over phase 
fluctuations, {e\p[iKZ]} = exp[-K^{Z}^ /2], and the fourth moment of the wave amplitude is 



r4(ri,r2) = 



1 



I 



d r\d r'2 exp 



/(r;-ri)-(r2-r^) 



-D^(r',)-D^(4)- 



D^{r[ + r^) ^ D^{r[ - r^) 



(27) 



Unfortunately, the form of eq. ( |27] ) in the strong scattering limit is not amenable to further simplification without explicitly specifying 
the form of the phase structure function. We thus note that the intensity variance is (A/^) = r4(0, 0) - 1 while the variance in the 
mutual coherence on a baseline r2 is (Ay^(r2)) = r4(0, r2) - exp[-D,/,(r2)]. 

The spatial power spectrum of the mutual coherence fluctuations across the observer's plane is obtained by making the replace- 
ment q = (r^ - V2)/r^ in eq. ( l27b . which allows the power spectrum to be identified as. 



WAviq; r2) = J 



d^r[ exp 



-iq ■ r'l - D^{r\ ) - D^(r2 + q^) + + 



■ r2 - qrj) 



2 2 

The corresponding power spectrum of the temporal intensity fluctuations is obtained by setting r2 = in eq. (l28T l: 



W/(w) 



2nv J 



d^r[dqsxp 



. DJr'. + qrl) DJr', - qrl) 



(28) 



(29) 



where the mass density power spectrum is assumed to be an isotropic function of q = (a>/v,q). Expressions for the intensity 
and mutual coherence fluctuations in the strong perturbation regime are only amenable to further simplification once the phase 
fluctuations have been specified; these are discussed in Appendix IB. 3 1 

Throughout the remainder of this paper we shall be mostly concerned with lensing in the limit of weak phase perturbations. It is 
shown in Sect.|7]that the large phase perturbation regime requires such large mass fluctuations as to be unlikely to apply even under 
the most extreme circumstances. 



4. Lensing by a power law mass spectrum 

In this section we examine the intensity fluctuations caused by a power spectrum of mass surface density fluctuations. This plausibly 
applies to lensing by small scale structure in cold dark matter if is clustered heirarchically on small scales. The power spectrum 
of dark matter fluctuations on very small scales is unknown, but it is not unreasonable to expect that it may follow a power law to 
sub-pc scales (e.g. Diemand, Moore & Stadel 2005). (The alternate possibility is that the matter is distributed in discrete clumps; 
lensing by such a distribution is examined in Sects. 5-6.) In the CDM paradigm the power spectrum of mass fluctuations follows 
a P(k) oc k"^^^ spectrum on the small spatial scales of interest here (Ch. 4 Padmanabhan 1993), where o- = 3/2 corresponds to an 
initially scale-invariant spectrum. The spectrum of fluctuations cuts off' at the free-streaming scale (Padmanabhan 1993; Peacock 
1999) which, for dark matter particles of rest mass energy E^, occurs at a spatial scale 

For cold (£■„, > 1 GeVc"^) dark matter particles the is well below the scales typically probed by gravitational scattering for wave 
frequencies v ^ 1 Hz. 
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Notwithstanding the distribution of dark matter on small scales, there are two additional reasons to consider lensing by a power 
law spectrum of mass fluctuations. The first is that it applies to lensing by any astrophysical turbulent plasma or gas. It is now well 
established that both the ionized interstellar medium of our Galaxy and the HI in our and nearby galaxies follow power law spectra 
over a large range of spatial scales which encompasses those of interest here, and there is every reason to expect this to apply to other 
systems (e.g. Armstrong, Rickett & Spangler 1995; Dickey et al. 2001; Stanimirovic & Lazarian 2001; Braun 1999). Although the 
lensing contribution from gas is likely to be small under most circumstances, it is nonetheless interesting to consider its contribution 
from a formal standpoint. It may also apply to lensing by a turbulent proto-galaxy. The second reason is that an examination of 
the effects of lensing from a power law mass spectrum provides a useful test-bed in which to develop an intuitive understanding of 
lensing in the long-wavelength regime applicable to gravitational radiation which can then be applied to more complicated lensing 
situations. We return to this point at the end of this section. 

We consider mass fluctuations with the following power law distribution, 



^ml' qo<q<q 



mm? 



f 2;(q) = QO i q-f^ ^min, <q< ^max, (31) 

( 0, q> q„ 



/max- 



This simple description is suflicient for all cases of interest here. Over the range of scales pertinent to lensing in the present case 
the mass fluctuations are expected to be described by a single power law index. One can, of course, prescribe a more complicated 
spectrum (e.g. by changing the behaviour at the inner and outer scales), but only at the expense of extra complexity. We see from the 
results below that such extra complexity is unwarranted in most cases. The spectrum flattens to a constant at an inner wavenumber 
9inin = i'o ' > where Lq is interpreted as the outer scale of the density fluctuations. For shallower power spectra, j3 <2, the outer scale 
of the power spectrum, g'^ax = 'o'' ni^st also be specified to ensure that the mass density variance is finite. A cutoff scale is also 
introduced to ensure the convergence of certain integrals evaluated below, but it will become apparent that this scale is physically 
unimportan{3. Expressed in terms of the mass surface density variance, the constant of proportionality Qq is. 



Bn " min ' ^ ' 



where ^min <^ ^max is assumed. For qo ^ ^min the structure function takes the following form (see Appendix lAlt, 

DAr) = 4nQoK^ \r^^ [log - i] - r^-l^jlllS- - ,4_^£™_ ^ (,6 4^/^\| o < /? < 4. 

^ I 4 [ ^Uo / ySj 23+/T(2 + f) 128(j8-2) V ' ^ 

(33) 

It is convenient to cast this in the form. 



, „ ^2 / „ \min[2+/?.4] 

D(r)= — - - , 0<^<4, r<q^, (34) 

\'"ditt7 \rQ 



where only the two most important terms are retained. This separation is useful because inspection of Z^(ri , r2) (cf. eq. (fTST i) shows 
that the terms make no contribution to the intensity fluctuations. The term proportional to dominates for r <K q^^, so the 
transverse separation between points on the lensing plane over which the rms phase changes by one radian is 



1 

'"difF 



?max(f=f'^, 0<,6<2 



, (35) 



However, since the term makes no contribution to intensity fluctuations, the salient quantity when considering intensity fluctua- 
tions is ro, which defines the length scale over which phase fluctuations that contribute to change by one radian rms. 



'"0 = 



25/4|/:|-i/2(Al2)-i74^ 2</3<4 



0</3<2 ^3^^ 



It is the magnitude of this quantity, and not rdift, which discriminates between weak and strong scintillation regimes. Strong scintil- 
lation corresponds to the regime ro <s: rp. Since it is far harder to satisfy the criterion ro « rp than rditf <K rp, the scattering can be 
often be classified as weak even when the relative phase delay between two points separated by a Fresnel scale on the lensing plane 
greatly exceeds unity. Thus the weak scattering condition holds over a larger range of scattering conditions than in other applications 
of scattering physics (see, e.g., Narayan 1992). 



^ Note that the form of the power spectrum differs slightly to that considered in M04, where the power spectrum of mass density (not mass 
surface density), Op, fluctuations was specified directly. The previous formulation is too simple for treating intensity fluctuations, as it included no 
power beyond the outer scale. 
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4.1. Intensity fluctuations 

The magnitude of intensity fluctuations in the regime of weak scattering is evaluated by inserting eq. dSTT l in eq. (fTST l. 



I dqq sin +2n \ dq q sin 



(37) 



In the Umit ^min'"F < 1 the modulation amplitude depends on the steepness of the mass power spectrum, with a change in the nature 
of the solution at y6 = 2. By recasting K^Qq in terms of ro, the scale length over which phase fluctuations responsible for focusing 
vary, we see that the magnitude of the intensity fluctuations in the regime of weak scattering depends only on the ratio rp/ro. 



25(^)\ 2<p<4. 



, 22+/^r(2 + f)sin(f)(^P, 0</?<2, 



Furthermore, the condition for weak scattering, rp < ro, implies that the intensity fluctuations in this regime are less than or 
comparable to unity. The amplitude of the intensity variations may be understood physically by regarding the phase delay imparted 
by the mass distribution as a perturbation to the phase curvature over the coherence region, of radius rp, that contributes to the 
intensity on the observer's plane. The rms intensity deviation is ^ feff//, where D^s - DlDls /Ds is regarded as the eff'ective 
distance of the observer behind the lensing plane and / is the focal length of the phase perturbation. The focal length of a lens that 
retards the wavefront by a phase A(p over a distance r from the lens centre is nk/Acp. The rms phase curvature across the Fresnel 
scale is D^ir^y^^, where we consider only the non-quadratic contribution to the phase structure function because the terms do 
not give rise to intensity fluctuations. The rms intensity deviation is thus (rp/ro)™'"'^'^''^^'^'. 

For very low frequency waves, v ~ 10"^ Hz, it is conceivable that the Fresnel scale becomes comparable to or even exceeds the 
outer scale of the power spectrum. For completeness we include the expression for the intensity variance in the limit ^min^p > 1, 

{Af}^n^K^Qoq-J^y,+0{q-^^^ (39) 

The intensity variance in this instance depends even more strongly on the outer scale of the power spectrum than in the q'min'"p < 1 
case. 

The foregoing discussion is related to the amplitude of the intensity fluctuations. However, several other observables are of 
interest, in particular the power spectrum of intensity fluctuations and the amplitude and power spectrum of fluctuations in the 
mutual coherence. These quantities are considered in detail in Appendix IbI 

Although the above results are derived using a specific model of the mass fluctuations, they lead to a physical understanding of 
the lensing which should apply more generally. The amplitude of the intensity fluctuations has a simple physical interpretation. It is 
governed by the phase curvature induced by the mass fluctuations across the coherence scale, which is set by the Fresnel length, and 
intensity fluctuations arise due to the focussing and defocussing of radiation across this scale. Accordingly, most of the fluctuation 
power occurs on timescales ^ rp/v, the time on which the line of sight traverses the Fresnel scale (see Appendix IB] for details). The 
utility of these results is demonstrated in the context of stellar microlensing below. 



5. The contribution from stars and other discrete objects 

5.1. Single lensing population 

Equations ( fTSl l and (l23T l demonstrate that the effect of lensing can be ascertained directly from the power spectrum of the mass 
fluctuations. Here we consider the effect of stellar microlensing by writing the power spectrum of the mass surface density in terms 
of the stellar mass profile and the spatial distribution of stars across the lensing plane. A simple prescription for determining the 
power spectrum for a collection of stars of identical masses was introduced in Macquart (2004), founded upon the work of Melrose 
(1996), but here we generalize the formalism to take into account the clustering properties of the stars and a distribution of stellar 
masses. Consider first the case of lensing by a collection of identical stars with positions r, on the lensing plane. The mass surface 
density takes the form, 

Nm 

2:(r) = ^/M(r-r,), (40) 

where /M(r) is the projected density mass of a star of total mass M centred on the origin. One is free to substitute any particular 
form for the stellar mass density profile, but for our purposes it suffices to approximate the stars as constant-density objects of radius 
R. The projected surface density associated with such objects is /M(r) - (3M/27tR^) r^, and the power spectrum of the mass 

surface density is 



sm(qR) - qR cos(qR) 



10 

The power spectrum of the surface density is constant for q < Vs/?"' 



M\l-^\+0{q^R'). (41) 
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The power spectrum of mass density fluctuations on the lensing plane which is related to the covariance of the mass surface 
density distribution, 

C2(r)=-i-2 r flf2qe'i>-(D^(q). (42) 

The covariance, C-^(r), of the mass distribution is computed by determining the Fourier transform of the mass surface density profile 
as follows: 

<S(r' + r)2(r')) = ^ J dh' Y(r' + r)S(r') = J d\ e''^■^£(q)|^ (43) 

where J[ is the area of the lensing plane. Since only fluctuations in the mass surface density across the lensing plane give rise to 
intensity fluctuations, one should strictly subtract the contribution of the mean density, (2), from the power spectrum. However, 
since the mean is obviously independent of position across the lensing plane it contributes only at the single point q = 0. Thus the 
subtraction of the mean density does not affect our results and the correction is subsequently ignored. 

The average power spectrum of the mass surface density is computed by separating the self (/ = j) contributions from the 
A^(A^ - 1) cross-term (/ + j) contributions and identifying cr - NjJi as the mean stellar surface density; 

1>2(q) = cr\fM{qt + cr(N - \)\fM{qt (e'l*''--^') . (44) 

The average over positions r, and appearing in the last term vanishes if the masses are distributed randomly over all space. 
However, this average is in general non-zero if the stars are clustered. Now if the distribution of stars is statistically homogeneous 
(i.e. it is wide-sense stationary), the joint distribution of object positions, /^(r,, r^), depends only on the difference, Ar = r, - ry, and 
is independent of the average position s - (r, -i- ry)/2. Writing /?(r,, r^) as a function of Ar only, the average over object positions is, 

(^"'"""■■'') = i / d'^vd^^,p{^v)e"^■^^ = (45) 

The function p{r) is regarded as the fraction of objects whose separation is between r and r + dr. Placing these averages back into 
equation (|45] l. the power spectrum of surface density fluctuations takes the simple form, 

<I>E(q) = cr\fM(q)\Hl+(rp2(q)] (46) 

This spectrum can be inserted directly into equations ( fTST i & (l23T l to determine the characteristics of intensity fluctuations due to 
lensing by an ensemble of identical stars. A specific model for stellar clustering is introduced below in which eq. ( |46] | is evaluated 
explicitly. 



5.2. Lensing by stars of different masses 

The generalization to stars of multiple masses is made by dividing up the stars into A^pop mass bins with the Nmj stars within each 
mass bin, j, located at positions r, j. The corresponding mass surface density takes the form 

Wpop ^M, 

= ('■-'■'■■>)■ (47) 

We proceed as above in computing the power spectrum of mass density fluctuations, which is now divided into three distinct sets of 
terms, (i) the self terms within each mass bin, (ii) cross terms within each mass bin (i.e. the contribution from stars of identical mass 
but different positions) and, (iii) cross terms from stars of different mass. Terms of the kind (i) and (ii) appeared in the single-mass 
case above, but the contribution from (iii) is new. The mass power spectrum is 

(48) 

In performing the average over star positions drawn from different mass bins we introduce the function pki{r), which measures the 
fraction of stars, of mass and Mi, separated by a displacement r. Thus the mass power spectrum takes the final form 

<l>2;(q) = Yj \fMM)\ [1 + o-MtPkk(q)] + o-M,.crM,/M,.(q)/M,(q)p«(q). (49) 

A- ktl 

Two properties of the model must be specified in order to compute the mass density power spectrum, (i) the mass distributions of 
the individual stars and (ii) the clustering properties. The effect of the latter is embodied in the probability distribution of interstellar 
separations, pki(r). Obviously if clustering is unimportant (i.e. the stars are distributed uniformly) this probability distribution is 
independent of r and all the cross-terms in eq. ( |49] l are then proportional [i.e. pkiiq) ^ S{q)] and are thus unimportant. The following 
subsection introduces a simple model to evaluate the effect of clustering on the mass power spectrum. The mass distribution is 
addressed subsequently. 
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5.3. Stellar clustering 

A statistical description of stellar clustering properties is important when the stellar surface density is high. Terms related to stellar 
clustering scale quadratically with the stellar surface density in the mass power spectrum, whereas uniformly distributed stars make 
a contribution that is only linearly proportional to the stellar density. Clustering properties enter into the power spectrum through 
the contribution of the cross terms o-M^iNuft - l)(exp[/q ■ (r,jt - ry.*)]) and o-MtiNM, - l)(exp[/q ■ (r/t - Tji)]} in eq. (l48T l above. 
The first cross term is due to the clustering between stars of identical mass, while the second represents clustering between stars of 
different mass. We make the simplifying assumption that stellar clustering properties depend only on stellar mass, so that stars of 
identical mass all possess identical clustering characteristics. 

The effect of clustering between stars of identical type is determined from the probability distribution of stellar separations 
between all pairs of stars on the lensing plane. Clustering information is embedded in the two-point correlation function, ^(r), in 
which the differential probability of finding two stars within areas dA i and dA2 respectively is 



dP = cr^[l + ^ix)]dAidA2, 



(50) 



where cr is the stellar surface density. It is unclear exactly how to parameterize stellar clustering, so we explore the following simple 
form which is found to apply to clustering between galaxies in cosmological surveys (Peebles 1980a): 



(51) 



where ro is regarded as the clustering amplitude to be determined either empirically or using simple physical arguments. Despite 
possessing the dimensions of length, it is not useful to interpret ro as a typical clustering length because the power-law form assumed 
for ^ for above is scale-free. Given the location of a star, the number of stars surrounding it within a radius R is. 



N{r<R)^o- d^r[l+^{r)] 



ttR^o- 



1 + 



■7 



(52) 



For < y < 2 the form of ^(r) in eq. ( fSTT l suffices to provide a physically acceptable description of stellar clustering over all possible 
separations, because the function N{R) is well-behaved over the entire range < ^ < oo. It satisfies the requirement that the excess 
probability of finding a star goes to zero at large separations provided y > 0. More importantly, it also satisfies the criterion that 
N{R) is finite at small R provided y < 2, which removes the necessity of imposing an arbitrary inner cutoff to hmit the value of ^(r) 
at small r. 

The number of stars within the thin annulus bounded by r + Ar and r is. 



dN . 
— Ar: 
or 



Inrcr 



1+1- 



Ar. 



(53) 



This quantity is directly proportional to the probability distribution of object separations. The probability distribution of object 
separations across the lensing plane is normalized by requiring that the total probability of finding the second object somewhere 
within is unity, yielding. 



p(r,e)drd0^ rdrde 



1 + 



ro 



nRi 



1 + 



Rn 



7 \ ro 



r<R^ 



(54) 



The first term in square brackets in the numerator represents the contribution for a uniform starfield, while the second term represents 
the additional contribution due to clustering. We use this distribution to compute the averages over object separations by assuming 
N » I and then taking the limit in which the lensing plane extends over a infinite area (i.e. so that Rout » ro). 



o-A?/e'q (r.-r,)\ ^ jjjj^ f dr f dee'"'"°"'p{r,0) 
\ ' Rout^o^Jo Jo 



r(i) 



< y < 2. 



(55) 
(56) 



Evaluation of the integral on the first hne above demonstrates that the contribution from the term representing a uniform stellar 
distribution is identically zero, so that only clustering contributes to this average over object positions. We also note that the average 
over object positions is positive since the ratio of the two F functions in square brackets is negative. 

Clustering between dissimilar populations is considered by generalizing the definition of the two point coiTelation function. The 
differential probability of finding a star of type 1 in the area dA\ and a star of type 2 in the area dAi is dP - 0-1O-2 [1 + ^i2(r)] dA^dAj. 
The two-point correlation function f i2(r) can be determined from the individual clustering properties of the two star types. Motivated 
by studies of cluster-galaxy clustering in cosmolog}|3, the two-point correlation function is approximated by the geometric mean of 



^ An example of clustering between dissimilar populations in cosmological contexts is the galaxy-cluster two point correlation, which can be 
compared against the galaxy-galaxy and cluster-cluster two-point correlation functions (see Peebles 1980b; Bahcall 1988). The galaxy-cluster 
two-point correlation function is empirically determined to contain two terms, the first of which is important at short scales and is due to the 
tendency of galaxies to cluster in clusters. In our treatment we assume that stars of a certain type have no tendency to cluster around stars of a 
different type, so this term is unimportant. The second term is related to the clustering of clusters themselves, and is important in our case. This 
contribution is given approximately by the geometric mean of the two-point correlation functions for objects of either type: ^ac ~ ^cc^gg- 
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the two-point correlation functions for objects of either type: ^/t/(r) = ^^t(r)'^^^;;(r)'^^ for clustering between stars of masses 
and Ml. With this generalization it is straightforward to calculate the contribution of the cross term as, 

crM,(NM, - 1) = -l^-yyn^^^cTM/J'^/J^lf-'-^^ (57) 

where ro,Mi and ro,M, are the clustering scale lengths of stars with masses Mk and Mi respectively. 

Although we have only explicitly considered stellar clustering in terms of a power-law two-point correlation function above, the 
generalization to other two-point correlation functions is obvious. 



6. Power spectrum of intensity fluctuations from stellar lensing 

Having obtained expressions for the power spectrum of mass density fluctuations across the lensing plane we proceed to evaluate 
the contribution of stellar microlensing on gravitational radiation. 



6.1. Lensing by a single population 

We first investigate the character of the intensity fluctuations due to only a single lensing population of stars. Inspection of eq. ( l44b 
shows that there are two distinct contributions to the intensity variance, the first from 'self-terms', which is proportional to the 
stellar number density, cr, and which we hereafter refer to as the 'incoherent' lensing contribution in light of the fact that the term 
is independent of the spatial distribution of the stars. The second term, whose contribution is proportional to cr^, is labelled the 
'coherent' lensing contribution in light of the fact that this term only contributes if the stars are clustered. These two contributions 
are most conveniently treated separately. When the stars are distributed randomly but uniformly over the lensing plane only the 
incoherent term contributes to the intensity variations. 

Intuitively one expects both the clustering properties and the Fresnel scale to strongly influence the nature of the intensity 
fluctuations. The lensing properties depend on the mass fluctuations contained within the coherence area determined by the Fresnel 
scale. The clustering amplitude determines the area over which the density fluctuation power spectrum scales quadratically with the 
stellar number density, rather than only linearly. The ratio of the clustering amplitude to the Fresnel scales thus bears on the character 
of the intensity fluctuations. For instance, we might expect that lensing by stars clustered on scales much smaller than the Fresnel 
scale would be equivalent to lensing by a uniformly distributed starfield, in which the power spectrum of intensity fluctuations scales 
only linearly with the stellar density. 



6.1 .1 . The contribution from uniformly distributed stars 



The Fresnel scale for the wavelengths typical of gravitational radiation is far larger than the radii of the individual lensing stars. 
Thus one expects the particular profiles of the lensing constituents to make a negligible difference to the magnitude of the intensity 
variance. For simplicity we employ the q < VS /? ' approximation to the power spectrum of stellar mass profile, |^(q)P - M^, 
and justify its use at the end of this section by examining the length scales of the fluctuations that contribute most strongly to the 
intensity variance. The intensity variance due to the incoherent term is. 



<A/^)ic = 4-o-M^K^ J cfq q'"' sin^ 



= n^o-M^K^rl. 



The power spectrum of temporal intensity fluctuations associated with the incoherent term is, 



WAcS)^ dql-^+qH 



1 - 



(58) 



(59) 



where we have now retained terms to second order in q in the spectrum of the stellar mass profile, |/(q)p. The temporal fluctuation 
spectrum is evaluated separately in the two regions q < r^'^ and q > r^K In the first region the sine-squared term in eq. (l23T l is 
linearised, and in the second its argument varies rapidly with q, so we approximate it by its mean value of 1/2. For co/v < rp' this 
approximation yields. 



InK^M^crrl 



•1, ,2 



1 



(60) 



The leading term inside the brackets dominates over the remaining terms by a factor larger than r^/R » 1. In the high-frequency, 
ll)/v > Tp', regime one has. 
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Fig. 2. An illustration of the importance of clustering in determining the contribution of stellar microlensing lensing. The top row 
depicts a starfield in which the probability distribution of stellar positions is uniform across the lensing plane, while the bottom row 
illustrates a distribution in which there is a contribution from clustering. Each left panel contains 10"* stars, and each right panel 10^ 
stars. As the stellar density increases the mass surface density begins to look relatively uniform for panels on the top row, leading to 
a uniform phase delay across the lensing plane compared to the case in which clustering is important. 



where again the leading term inside the brackets dominates. One thus sees that the temporal power spectrum is constant for w/w ^ 
rp' and then decreases as af^ at larger angular frequencies. Since the contribution to the intensity variance scales as (DWitS), 
structure with wavenumbers q ~ r^^ dominates the contribution to the intensity fluctuations, and intensity fluctuations occur most 
prominently on a timescale r ~ r-pjv. 

For the wavelengths typical of gravitational radiation the Fresnel scale far exceeds the radii of the individual lensing stars, 
this result also justifies our original approach in approximating the stellar mass profile |/M(q)P to lowest order in q only. This 
approximation is also valid for the lensing of radio-wave electromagnetic radiation at cosmological distances. However, at shorter 
wavelengths the approximations used above fail, where the radii of the stars themselves exceed the Fresnel scale, the above integrals 
need to be re-evaluated using a precise expression for the stellar mass profile. 



6.1.2. The additional contribution from clustering 

The foregoing results provide a lower bound to the contribution from stellar microlensing. However, the contribution can be much 
greater than if the stellar density in the lensing region is sufficiently large and the stars exhibit a tendency to cluster. Figure |2] 
demonstrates why one intuitively expects the contribution of clustering to dominate once the stellar density becomes large. For 
an unclustered stellar distribution the rms deviation in mass density increases only as the square-root of the density, whereas the 
contribution from clustered stars increases linearly. 

The contribution to the intensity variance from clustering terms is evaluated using the average over object positions in eq. ( |56] l. 
and again using the approximation |/(q)P - M^, 



(A/^)coh = SnK^ J qdqq "^sin^ 



2 2 \ 



o-A^(exp[!q- (rz-r^)]) 
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The importance of clustering in determining the lensing characteristics is illustrated by the fact that this contribution from the cross 
terms is a factor ~ crr^ (fo / r^Y larger than that from the self-terms. This result is easily understood when one considers that the extra 
number of stars within the Fresnel radius when clustering is important is N{r < rp) ~ nr^cr{r-p I rQy~^ . The stellar number density 
need not be very high for clustering to dominate the amplitude of the intensity fluctuations; one requires a stellar number density 
that exceeds one star over an area r^irolrpY . This area can be of order many square parsecs if the lensing occurs at cosmological 
distances. 

The temporal power spectrum is again evaluated by approximating the sine-squared function as q^r^/A- for q < rp' and as 1/2 
for q > r„^. The calculation of the spectrum is detailed in AppendixO and we merely list the results here: 



(\ f —4 — ( — 

r(-r)xi , 'iy > \ " ■ 
2 I 1 Q (y-3)(y-i) 1 U>Y ^ y-i 

[ °(y-2)(y-4) v\v) ' 'f 



= cr^M^K^r^y sin 1 1^ ) r(-y) X <{ " ^ ^ ''^^ ^ ,^ ^' ^ (63) 



The power spectrum scales as w'' in the low frequency regime but decreases sharply in the regime lo > vr^ . Most of the 

.-1 



fluctuation power, which scales as ojWioj), is therefore sharply concentrated on angular frequencies to ~ vr.^ 



6.2. Generalization to a distribution of stellar masses 

The generalization to a population of stars requires a specification of both the mass distribution and clustering properties. The mass 
distribution of stars is taken to be a power law with index so that if the total stellar number density is cr the surface density of 
stars of masses between M and M + dM is. 



dM ^ M 
— = Co- — 



dM 

Mo' 



(64) 



between lower and upper cutoffs M^in and M^ax respectively. The constant of normalization, C - (1 - j6)[(M,nax/A^G)' - 
(Mtnin/Mo)'"^]"', is chosen so that integration of ctm over the mass distribution equals the total stellar number density, cr. A 
practical specific choice is a Salpeter initial mass function (Salpeter 1955), in which /3 - 2.35, the mass distribution ranges between 
0. 1 and 125 Mg and the constant of normalization is C = 0.0603. 

If there are sufficiently many stars one converts to the continuum limit in which the sums over object masses in eq. ( |49] l are 
written as integrals, so the mass power spectrum becomes. 



3'5;(q) 



Co- 
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where ro is parameterized in terms of stellar mass. The radii of the lensing stars also depend on the mass, but this only changes 
the mass power spectrum at wavenumbers q ~ R^. This wavenumber remains sufficiently large even for the largest stars that it 
does not appreciably alter the lensing properties because, as shown above, the dominant contribution to the intensity variance arises 
from mass fluctuations on length scales ~ rp » R for the wavelengths typical of gravitation radiation. However, for lensing of 
electromagnetic radiation at wavelengths shorter than ~ 0.01 m this simplifying assumption fails, and the density profiles of the 
stars must be taken into account; this case is not investigated further here. 



6.2.1 . A specific model for tine mass and clustering distributions 

We introduce here a simple analytically tractable model which, though possibly unrealistic as a complete description of stellar 
clustering, is nevertheless useful as a context in which to explore certain properties of clustering germane to an understanding of the 
underlying physics of stellar microlensing. 

Our prescription for the clustering length scale is founded on the basis that long-lived stars diffuse far from their birthplace over 
their lifetime and thus cluster weakly, but higher mass stars are more short-lived and are thus highly clustered because they have 
less time to diffuse from their birthplaces, in which they are originally highly clustered. We thus postulate a clustering amplitude of 
the form. 



I MY 
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where the positivity of a reflects the fact that less massive stars are more uniformly distributed and thus possess smaller clustering 
amplitudes. The intensity variance is calculated by replacing the integrals over object masses in expressions (|58l ) and ( |62] i to give. 
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where the q < yfSR ' approximation to |/M(q)P is again employed. 

Similarly, to obtain the power spectrum of intensity fluctuations one makes the replacement, 
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The clustering contribution from stars of various masses depends on how the clustering depends on mass. It might be supposed 
that massive stars with short lifetimes diffuse relatively small distances from their natal molecular clouds, and thus cluster strongly 
compared to dwarf stars, which may diffuse long distances over their much longer lifetimes. For a distribution following a Salpeter 
mass distribution low mass stars make the dominant contribution to the intensity variance when a < 0.7, both low- and high-mass 
stars are important in the range 0.7 < a < 1.7, while high-mass stars make the dominant contribution to the intensity fluctuations 
when a > 1.7. 



6.3. Summary of results 

We provide here a brief summary of the important results pertaining to the lensing of gravitational radiation by stars derived in the 
preceding two sections. The amplitude of the intensity fluctuations is determined by the mass density power spectrum. For lensing 
by discrete objects, this power spectrum contains two sets of terms: (i) a set of self-terms which is linearly proportional to the 
object density, cr, and is independent of how the objects are distributed across the lensing plane and, (ii) a set of cross-terms which 
increases quadratically with the stellar density but is only important if the objects are clustered. The latter terms are important for 
gravitational radiation because the stars may cluster on length scales comparable to the coherence radius rp, so that a collection of 
clustered stars acts like a single lensing macrostar of mass MN. The lensing contribution from the same number of uniformly 
distributed objects is equivalent to an object of mass M y/N only. In a simple case in which stars are clustered across the lensing plane 
according to a power-law two-point correlation function ^(r) = (r/roy', the contribution to the intensity variance from clustered 
stars exceeds the uniform component by a factor (Tr^(ro/rpy. 

The effect of lensing by objects of different masses is considered by investigating the effect of lensing from objects whose mass 
distribution follows a power law. For uniformly distributed stars the intensity variance is dominated by the high mass cutoff of the 
object distribution provided that the mass distribution falls less steeply than M^^. Thus lensing by a starfield that follows a Salpeter 
initial mass function (i.e. a distribution proportional to ~ M"^ 35-) jg dominated by the high mass cutoff when clustering is negligible. 

The power spectrum of temporal intensity fluctuations is also evaluated for both uniformly-distributed and clustered objects. It 
is shown that in both cases the largest contribution to the fluctuations occurs on the timescale required for the objects on the lensing 
plane to drift a distance ~ rp across the line of sight. The lensing is also most sensitive to mass fluctuations on this length scale. 



7. Signatures of lensing and magnitude estimates 

We present here several simple scaling relations that allow the magnitude of intensity fluctuations to be computed for any given 
lensing geometry and mass fluctuation amplitude. We defer to a later a work a more rigorous analysis that takes into account the 
distance ranges out to which the various types of source might be observed by current and future detectors. 

The magnitude of the intensity variations depends on the strength of the scattering, which is set by the amplitude of the mass 
power spectrum relative to the Fresnel scale. It is convenient to define an effective distance to the scattering material, Dgff = 
DlDls IDs so that the Fresnel scale takes the form rp - [AD^ffjil + Zl)]'^^. When considering scattering by an intervening system 
this effective lensing distance is of order the angular diameter distance to the system itself. In terms of values normalised for values 
typical of scattering by matter from a source at a cosmological distance (e.g. a galaxy), the Fresnel scale for lensing at a distance 

£>efF is. 
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7.1. Lensing by a power law spectrum of mass fluctuations 

For a power law spectrum of mass fluctuations the scale length of the phase changes that give rise to intensity variations is 

\-l/3 



'"0 = 



665(l.z.)-/3(^r(^r(^)-'pc,/..l, 
213(l.z.r'/^(^r^^(^)'%c, 2</3<4 



(71) 



where, in order to ensure tiiat the numerical constants do not themselves depend on /3, the case = 1 is chosen to represent the 
regime < /3 <2. Geometric phase delays thus dominate phase delays induced by inhomogeneities in the mass distribution, so that 
the intensity fluctuations occur in the weak perturbation regime. 

In calculating the magnitude of lensing effects on gravitational radiation it is important to recognise that gravitational wave 
detectors measure the wave strain directly, so it is most appropriate to consider fluctuations in the magnitude of 0, where even a 
small intensity variance implies potentially measurable fluctuations in the gravitational wave amplitude. We express the variability 
in terms of the rms wave amplitude fluctuation, {6\<p\^y^^. In the regime of weak fluctuations it can be shown using the Rytov 
method that the intensity fluctuations are log-normally distributed (Tatarski 1967; Fante 1968) so that here, where the fluctuations 
are very small, this is well-approximated by a normal distribution (see Goodman (1985) and references therein). The rms wave 
amplitude fluctuation, {6\(p\^y^^, is thus expressed in terms of the intensity variance as (2/;7r)'^'*(A/^)'^'*. We also explicitly write 
observational quantities in this section in terms of the mean wave strain of the source, which is the value that would be 
measured in the absence of lensing effects. (In previous sections this normalization was made implicit by assuming the source to be 
of unit amplitude (cf. Sect. 2).) The rms wave amplitude fluctuation is thus 



4x1/4 



= 0.012 I —^^^^^^ I I -p-^ I , 2 < y3 < 4. (73) 

\iM2pc-4; \iGpc; 

If the wave amplitude fluctuations diff'er from a Gaussian distribution the results will diff'er only slightly from those in eqs. ( l72t & 
( f73T l. This is because the rms of the wave amplitude fluctuations differs by only a value of order unity to the power of one quarter 
if the phase fluctuations deviate from a normal distribution. We note that the estimates in eqns. ( |72] | & ( l73T l are similar to those 
computed by Takahashi (2006), who considered fluctuations in the magnitude and phase of (p directly. 

The magnitude of lensing-induced variations depends on the mass surface density variance. This quantity may be very large 
if the line of sight intersects a massive structure, such as a galaxy cluster. We therefore compute a worst-case upper limit to this 
quantity based on measurements of the mass surface density in galaxy cluster environments using weak lensing measurements. A 
suitable fiducial mass density E a; 140 pc"^ is derived from weak lensing measurements of the cluster CI 1 358-1-62 (Hoekstra et 
al. 1998). If we assume that the bulk of this matter is distributed inhomogeneously over the extent of the cluster one has (AZ^) ~ 
2.0 X lO'* M| pc The rms wave amplitude fluctuation is thus a; 15 % for the fiducial numbers used in the scaling relations ( l72b & 
( |73] ) above. 



7.2. Lensing by stars 

An estimate of the stellar density on the lensing plane is crucial in estimating the magnitude of the stellar lensing effects. The 
density also determines the number of stars that contribute to the lensing at any one instant. The number of stars that contribute 
simultaneously to the lensing is determined by the coherence area, whose radius is comparable to the Fresnel scale, which depends 
on the effective lensing distance. For lensing within our Galaxy one has Dgir ~ AL/2 ~ 5 kpc for a medium of thickness AL ~ 10 kpc, 
whereas for lensing of an object at cosmological distances by an intervening system, the fiducial effective lensing distance is of order 
1 Gpc. For lensing within the Galaxy the number of stars that contributes to the lensing is, 

A^MW~0.2v-'(-^l(-^f , (74) 



Ipc ■'/y 10 kpc 

for a stellar volume density . The number of stars that might contribute to lensing by an intervening system for a source located 
at cosmological distances is much greater, 

«.,c^.«'o...-.-(^)(^)(^). 

Thus the statistical approach adopted here is nearly always applicable to intergalactic scattering, but is only useful for Galactic 
lensing at frequencies v ^ 10"^ Hz. The following discussion is restricted to regimes in which there are sufficiently many lensing 
objects that a statistical approach is useful. The effect of lensing by only a few discrete objects has been considered extensively 
elsewhere (e.g. Chang & Refsdal 1979, 1984; Kayser, Refsdal & Stabell 1986; Schneider & Weiss 1987; Wambsganss 1990; Lewis 
et al. 1993). 
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The simplest feasible lensing scenario consists of lensing by a large number of identical stars. This approach presents a service- 
able description of lensing by dwarf stars whose spatial distribution is likely to be devoid of clustering because their long lifetimes 
allow them to diffuse homogeneously throughout a galaxy. The contribution from such stars is described by eq. dSST l. which gives. 



0.019(1 +zz,y^^(—) I I {tttt^] . (76) 



<|<^|) ■ ^ ' llHz/ \lGpc/ \ 100 stars pc-2/ \o.lM, 

The same lensing contribution is obtained if the lensing plane were instead populated by 1 stars with surface density 1 starpc"^. 
The generalization of this result to lensing by a distribution of stars of various masses shows that the low mass end of the distribution 
dominates only if the index of the power-law distribution falls more steeply than M^^. Most expected mass distributions are shal- 
lower than this, so the dominant contribution to the lensing amplitude arises from the upper cutoff in the stellar mass distribution. 
Taking the distribution to follow a Salpeter initial mass function with O.IM© < M < 125M0 and /3 - 2.35, eq. ( |67] ) shows the rms 
amplitude fluctuation is, 

(101) llHz/ \lGpc/ \ 100 stars pc-2/ 

When the stars are clustered there is an additional contribution to the intensity variance. The magnitude of this contribution 
depends on the clustering amplitude, which is poorly constrained a priori. It is convenient to parameterise the clustering amplitude 
in terms of an over-density parameter, x, which quantizes the factor by which the stellar density exceeds the mean (uniform density) 
over some fixed separation. In this way, the number of objects within a region of radius X is N{r < X) = (I + x)o'kX^. We choose 
a fiducial distance X = 1 pc and accordingly label the overdensity parameter ;^'ipc. The additional clustering contribution to the 
intensity variance for a single clustered lensing population is thus, from eq. ( |62] i, 

} UHz/ UGpc/ \ 100 stars pc-2j \0.1Mo/ I 1 / ' 



(78) 

where the numerical constant reflects the choice y - I (cf. eq. (ISTTl). The rms wave amplitude for other values of y is obtained by 
noting that the amplitude varies slowly, between 0.107 and 0.128 for y between 1/2 and 3/2 for the fiducial parameters used in 
eq. ( l78l l. The generalization of this estimate to a distribution of masses requires assumptions about the strength of clustering as a 
function of mass. Taking a = 1 as the index which characterizes the dependence of the clustering properties on mass (cf. eq. (|66]|). 
it remains only to determine the clustering amplitude for a stars of a single mass. If the number overdensity is;^'ipc for stars of mass 
Mc, eq. ( |67] ) shows that the rms wave amplitude is, 

{\cf>\} ■ ^ ^' llHz/ \lGpc/ \ 100 stars pc-2j I 1 I [lOMsj ' 

(79) 



where y = 1 is assumed. 

In the foregoing discussion the contributions from uniformly distributed and clustered stars are treated separately. The contri- 
butions from both are summed to deduce the total intensity variance. Similarly, the total rms fluctuation in the wave amplitude is 
determined by summing the incoherent and coherent parts raised to the fourth power (i.e. (6\<p\^}^ - ((^101^)?/^)'* + ((<5|0p)'^^)"^). The 
total rms is thus well approximated by the larger of the two individual contributions. 

Some remarks on the applicability of our formalism to the lensing of electromagnetic radiation are in order. To lowest order in 
the wave amplitude the propagation of both electromagnetic and gravitational radiation are subject to the same wave equation. The 
main distinguishing feature between electromagnetic and gravitational radiation is the wavelength and hence the magnitude of the 
Fresnel scale. Thus the coherence area and hence the number of stars that contribute to the lensing simultaneously is smaller for 
electromagnetic radiation. The smaller value of the Fresnel scale renders the effects of stellar clustering much less important for 
electromagnetic radiation. The contribution to the intensity variance from clustering scales as /l'"''^^. 

However, the presence of fewer stars within the coherence area does not invalidate the statistical approach adopted here. This 
is because the ensemble-average quantities apply equally whether one averages over space or time. Obviously a statistical treat- 
ment is unnecessary if only a single star contributes to the lensing at a time, but the ensemble average values calculated here are 
nonetheless correct provided only that the average is of sufficient duration to encompass many lensing objects drifting past the 
line of sight. Although the formalism developed here is applicable to electromagnetic lensing, we caution that the final results 
for intensity variances and power spectra make the approximation that the Fresnel scale greatly exceeds the radii of the lensing 
objects and thus that the shapes of the individual lensing objects are unimportant. This is obviously an excellent approximation 
for gravitational radiation. However, this approximation fails at optical wavelengths for lensing at cosmological distances. More 
generally, the mass profiles of the individual stars must be taken into account for lensing at an effective distance Deif at wavelengths 
A ^ lOQ (Deft/l Gpc) nm. When this assumption fails the expressions calculated here for the intensity variance and temporal power 
spectrum must be recomputed, taking into account the particular mass surface density profiles of the lensing objects (cf. eq. (14111). 
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7.3. Signatures and timescales of tensing 

Although there are several signatures which could potentially identify radiation that has been lensed by small scale structure, several 
mitigating factors make their unambiguous detection problematic. The most obvious characteristics are that the fluctuations are 
stochastic in nature and that the power spectrum of the variations contains a number of prominent features related to the wavelength 
and lensing geometry. The temporal power spectrum exhibits breaks at several characteristic frequencies, most notably at w = Wp', 
beyond which it steepens. Analysis of the lensing by both stars and a power law distribution of material shows that the intensity 
variance is dominated by fluctuations on a timescale r x rp/v: 



It is difficult to detect temporal fluctuations because most of the relevant fluctuation timescales are too long to be detected. For a 
source or lens traversing the line of sight at even ~ lO^kms"', the expected fluctuation timescale is > lO-'yr for the estimate of 
the Fresnel scale in eq. ( fTOl i at a frequency of 1 Hz. Shorter timescales are, of course, possible if the source moves at relativistic 
speeds, but this is unlikely to occur in practice. The implication here is that any source lensed by small scale structure will not vary 
appreciably due to lensing, despite the fact that its intensity is strongly altered by lensing. This presents a problem for gravitational 
wave detectors hoping to remove the effects of stochastic lensing, because the detectors will not be able to average over multiple 
fluctuations and recover the intrinsic source intensity. 

The frequency dependence of the intensity fluctuations depends on the slope of the underlying mass density power spectrum. For 
steep, /? > 2, power law spectra the mean intensity variance is independent of frequency, while for shallower spectra the intensity 
variance exhibits a moderate, v'"^^^, dependence. The frequency dependence of the intensity variance for lensing by stars is in the 
range v'^^ to v^^''. Any frequency dependence is potentially observable in chirped (e.g. inspiralling) gravitational wave sources. 
The amplitude of the stochastic intensity fluctuations would decrease with time, as the emission frequency increases. However, 
identification of the frequency dependence due to lensing is complicated by the fact that the scaling only applies to the ensemble- 
average intensity variance, whose estimation requires intensity measurements over several fluctuation timescales at each frequency. 
This appears to render practical measurements of the frequency dependence unattainable. 

8. Conclusions 

A formalism for treating lensing by (i) objects distributed according to a power law spectrum, and (ii) stars and other discrete objects 
is considered in the long wavelength limit applicable to gravitational radiation. For the large wavelengths typical of gravitational 
radiation and the large distances characteristic of lensing at cosmological distances the coherence area, of radius ~ rp, over which 
objects influence the observed wave properties exceeds many square parsecs. Since this area is likely to contain a large number of 
objects for any line of sight which intersects a galaxy on its way to Earth, a statistical treatment is necessary to compute the effect 
of microlensing. 

For a given amplitude of mass density fluctuations the amplitude of the wave strain increases with the effective lensing distance 

1 /2 

as ~ D^^ . Thus the effect is most pronounced for sources at cosmological distances in which the line of sight intersects a prominent 
collection of mass, such as a galaxy or proto-galaxy. The effect is thus much smaller for the lensing of nearby, Dgfr ^ 100 Mpc, 

1 /2 

sources for two reasons. Namely, the Fresnel scale (for which rp oc DJ^. ) is small, and the radiation is less likely to encounter 
an intervening galaxy. In the absence of lensing by intervening systems only material within the source's own galaxy is likely to 
contribute to the lensing, for which the effective lensing distance is only of order the size of the galaxy (i.e. tens of kpc at most). 

For lensing at cosmological distances stellar microlensing induces rms wave strain fluctuations of order 1 - 2% for unclustered 
stars and 5 - 10% for clustered stars with stellar densities on the lensing plane of ~ lOOpc"^. However, it may be expected that dark 
matter, which may also be highly clustered, and whose contribution to the total mass of a galaxy exceeds that of luminous matter, 
may make a much stronger contribution to the intensity fluctuations. Lensing by intervening small-scale fluctuations in the dark 
matter can induce 2-15% rms variations in the amplitude of gravitational radiation. 

Temporal intensity variations caused by movement of the lensing material across the source-observer line of sight occur on too 
long a timescale to be observable. Most fluctuation power is associated with mass variations on scales comparable to the Fresnel 
scale, and most of the fluctuation power is thus on timescales exceeding rp/v, for a transverse speed v plausibly at most lO-'km s"' . 
The Fresnel crossing time for a lensing speed of 10^ km s~' exceeds 1000 v ''^years. 

Intensity variations due to lensing provide the opportunity of exploring the nature of mass fluctuations in lensing systems over a 
large range of scales. Lensing of ~ 1 Hz gravitational radiation at cosmological distances is sensitive to mass fluctuations on scales 
of order a parsec. Radiation detected by LIGO at frequencies ~ 10 kHz is most sensitive to structure on 0.01 pc scales, and radiation 
detected at the bottom end of the LISA band, v ~ 10"^ Hz is sensitive to structure on ~ 300 pc scales. The fluctuation power is 
attenuated if the lensing medium contains no power on these scales. Thus investigation of the gravitation wave lensing constitutes a 
sensitive probe of the dark matter power spectrum over a range of relatively small scales. 

We have also considered variations in the mutual coherence of the wave strain, the quantity (0(r')0*(r' + r)), which is of interest 
because it is the most basic estimator of the effect of gravitational lensing on the wave strain. The degree to which the mutual 
coherence is a useful statistic is estimated by calculating the variance of this quantity. We find that the amplitude of fluctuations in 
the mutual coherence is comparable to that expected for intensity fluctuations except when the mutual coherence is measured on 
baselines r large compared to ro, the scale over which the fluctuations in the phase curvature occur. This scale is very large so that, 
in practice, the mutual coherence fluctuates in the same manner as the intensity, (p(p*. 
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Appendix A: The Structure Function for a Power Law Spectrum 

The phase structure function corresponding to this spectrum is 
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We expand the hypergeometric functions for q^c^^r « 1, using 
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Appendix B: Characteristics of lensing by a power law mass spectrum 

B. 1 . Power spectrum of intensity fluctuations 

The power spectrum of temporal intensity fluctuations provides information on the relative contributions of various scale lengths to 
the intensity variance. The temporal power spectrum is given by, 
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There are three regimes over which the behaviour of the power spectrum changes character for the case q„^i„ < : (i) co/v < q'min < 
rp', (ii) q^in < a)/v < r^^ and (iii) a>/v > rp'. If instead Fresnel scale exceeds the outer scale (^min > '"p') then only cases (ii), 
oj/v < ?"p', and (iii) are applicable. 

The limiting behaviour in case (i) is evaluated by linearising the sine functions appearing in the integrals. For cases (ii) and (iii) 
only the second integral in eq. dB.lt contributes to the power spectrum. For case (ii) we again linearise the sine function, while in 
case (iii) the argument of the sine function is much larger than unity and we approximate it by its average value of 1 /2. This gives 
the following asymptotic behaviour for /3 ^ 3/2: 
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Thus, for a given Fresnel scale, we see that the power spectrum consists of a flat portion at low frequencies, a portion over which the 
spectrum decreases as w'"^ at intermediate frequencies, and a portion that increases even more steeply, as oj^^^^, at high frequencies. 
In Figure 2 we plot both a numerical integration of eq. ( IB. II ) and its asymptotic behaviour as derived in eq. ( IB.2l i. 

The part of the power spectrum that corresponds to short timescales, a)/v > rp', contains an oscillatory component that varies 
approximately as sin^ o/r^/v^ and which is not modelled by the asymptotic solution. The oscillations become increasingly rapid at 
higher frequencies. These oscillations are an inherently wave-optics effect attributable to constructive and destructive interference 
caused by oscillations in the degree of coherence due to the geometric phase delay term (i.e. the first term in the exponential in 
eq. (|2]i). They may be understood by recalling from elementary optics that scattering from a plane in which the lensing-induced 
fluctuations are small compared to the geometric phase delay shows a number of concentric bright and dark rings. This situation is 
analogous to the familiar pattern displayed in Newton's rings, in which a convex interface placed on a flat surface and illuminated 
with monochromatic light displays a bright central image suiTounded by a number of concentric rings. The dark and bright rings 
are due respectively to destructive and constructive interference between radiation reflected from the upper and lower surfaces. 
The central bright peak is due to structure within the first Fresnel zone, the second peak is related to the second ring of coherence 
surrounding the phase centre, and so on. Thus we see that the first bump observed in the power spectrum at oj/v ~ r^^ is due to 
interference within the first Fresnel zone of the scattered image. The next peak in the power spectrum is due to interference caused 
by the second Fresnel zone, and so on. The period of the Fresnel oscillations in the power spectral domain depends only on the 
Fresnel scale and the effective lensing velocity. 

The most important feature of the intensity power spectrum, however, is that most of the fluctuation power occurs on long 
timescales. The contribution to the intensity variance over an interval Aco scales as ~ <jjWj(a>)Aa). Thus the peak contribution for a 
mass spectrum steeper than /3 - 2 occurs atco ~ vq„^i„. The contribution over the range a>/v - [0, <7min] exceeds that over the interval 
[^min, '■p'] - This implies that an observation over a time interval t > Lq/v is required to obtain an average source intensity close to its 
ensemble-average (intrinsic) intensity. If the lensing-induced intensity variance is large, this means that no practical measurement 
of the average intensity will approach its ensemble average value for the mass spectrum outer scales and effective lensing speeds 
likely to be encountered in reality. The situation is less dire for shallower,/? < 2, spectra, where the peak contribution to the intensity 
fluctuations occurs on shorter timescales r ~ rp/v. The divergence of observational quantities from their ensemble-average values 
is explored further in the discussion below, since it bears ramifications for the suitability of gravitational wave sources as standard 
sirens. 



B.2. Fluctuations in the mutual coherence 



The variance of the mutual coherence is evaluated assuming r2 < rp, which allows us to expand the sine-squared function in the 
integral of eq. ( |25] l to second order in q ■ r2. The variance in the visibility measured on a baseline r2 is. 
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Fig.B.l. An illustration of the behaviour of the power spectrum of temporal intensity fluctuations, Wiioj), for lensing through 
a power law spectrum of mass fluctuations. The solid line is obtained by a numerical integration of eq. ( IB. lb with inner scale 
?min - 0.01, Fresnel scale rp - 1 and index [3-3. The dashed line is the asymptotic solution of eq. ( IB. 21 ). 



The first term is recognised as the intensity variance, while the second term appearing inside the integral is identically zero for any 
isotropic mass power spectrum. For the power law mass distribution considered above the third term is evaluated and expanded to 
lowest order in to give 



(B.4) 



We calculate the power spectrum of fluctuations in the mutual coherence to estimate the timescale on which it fluctuates, and 
ascertain whether a measurement of the mutual coherence is likely to correspond to its ensemble average value. The investigation is 
restricted to the case r2 || v since this is the most relevant to observations with a single detector. The detector measures V{r2 = vf) 
by measuring the scattered wavefield as it drifts past the detector Thus one is most interested in fluctuations in the mutual coherence 
when the baseline over which it is measured is also oriented parallel to the relative lensing velocity. The power spectrum of mutual 
coherence fluctuations is evaluated in the same manner as the intensity power spectrum to find. 
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Once again, we see that most of the fluctuation power is on the largest scales. This implies that an observation over a time interval 
T > Lq/v is required to obtain an average visibility close to its ensemble average value for /3 > 2 spectra, while integration times 
exceeding r ~ rp/u are required when </3 <2. 



B.3. Strong phase perturbations 

Fluctuations in the strong scintillation regime are evaluated using the phase structure function, eq. ( l33T l, rather than the power 
spectrum directly. Strong scattering only applies when KZ(ri,r2) > 1 over the effective region of integration in ri and r2 in eq. 
ST3[ . The quadratic part of the structure function makes no contribution to Z and thus to the intensity fluctuations. Physically this 
is because the quadratic term is associated with (linear) phase gradients which possess no curvature. Since only the next most 
important term contributes to intensity fluctuations, this regime only applies to scattering in which the mass density variance is 
exceptionally high or the Fresnel scale is large. Although it is demonstrated in the next section that strong scintillation is unlikely to 
occur even under the most extreme circumstances, the behaviour of the intensity fluctuation is still of formal interest. The analysis 
is divided into two domains: (i) < y6 < 2, where the part of the phase structure function that contributes to intensity fluctuations 
scales as r^^^, and is a factor (r/Lo'f smaller than the quadratic term, and (ii) 2 < /? < 4 where the dominant term scales as r* and is 
a factor of (r/Lo)^ smaller than the quadratic term. 



B.3.1./?<2 

A number of authors have investigated this case in the context of scintillation (Jakeman & Jefiferson 1984, Rino & Owen 1984, 
Goodman & Narayan 1985), so we merely summarize these results and recast them in the terms of the notation employed here. In 
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the regime ro <K rp applicable to the strong perturbation case the intensity variance is independent of the scattering strength and is 
asymptotically equal to 



2#TT 
2-P 



1, 



0.2 ^yS^ 1.8, 



0.2 g/? ^ 1.8, 



(B.6) 
(B.7) 

(B.8) 



where the subscripts 'diff ' and 'ref denote the character of the intensity variations as diffractive and refractive respectively. In 
practice these results are valid for rp/ro ^ 5 (Jakeman & Jefferson 1984), but the asymptotic solution is unreliable for values of /3 
near both and 2. 

Refractive variations occur on the timescale. 
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while diffractive variations occur on the much shorter timescale. 
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Refractive variations are caused by phase curvature over the entire image of the scattered source as projected on the lensing 
plane. There is a large amount of power on large scales because the phase power spectrum scales sharply as q^'^^^. Following 
the argument advanced by Goodman & Narayan (1985), the root mean square phase curvature per logarithmic interval scales as 
~ QoK^q^^d(}ogq), and structures on length scales q focus at a distance z ~ k/[QoK^q^^^d{logq)y^^. Thus the largest fluctu- 
ations which focus on the observer's plane, which are also most important due to the steep phase spectrum, are those at spatial 



wavenumbers ^,-ef ~ r'^^^ 2)^(^2+/s)(2 /3)^ ji^i^ gives rise to lensing on timescales Tref 
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Diffractive variations are caused by interference between radiation arriving from different parts of the lensed image. This pro- 
vides a physical argument for the timescale of diffractive changes. The interference pattern changes on a timescale over which the 
relative phase delay between two patches separated a distance ~ q^^^ on the lensing plane changes by order one radian. If both 
patches move a distance Ax across the lensing plane the phase delay changes by ~ + Ax)^ - (Ax)^]/2rp. Thus variations occur 
on a timescale t ~ r^q^gf/v, which explains the magnitude of Tdiff above. 

The variance of the mutual coherence is computed by Goodman & Narayan (1989) [cf. their eqs. (4.2.2) & (2.5. 12)] who find. 
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where K is a constant that is not of interest here. The variance in the mutual coherence decays over a much larger scale than the 
mean visibility itself. The mean visibility decays exponentially quickly on a scale length rdiff. Recall that, unlike for intensity or 
visibility scintillations, the quadratic term in the structure function makes the dominant contribution in determining the ensemble 
average visibility (i.e. {V{r)) = exp[-0.5(r/rdifr)^]). However, the mean square visibility varies on a much larger scale set by 
ro » /"diff ■ Simulations of scattering from such steep spectra (Narayan & Goodman 1989) attribute this difference to the fact that the 
instantaneous lensed image has a fractal patchiness, meaning that an instantaneous image of the source appears fragmented into a 
number of highly-modulated subimages. The ensemble-average image is necessarily smooth because it averages out all the effects 
of the large (linear) phase gradients and thus image wander, which are the primary cause of this fragmentation. 



B.3.2. 2 < yS < 4 

For yet steeper mass density power spectra the part of the phase structure function that contributes to intensity fluctuations scales as 
r'*, and the dominant contribution to the fourth-order moment of the wavefield reduces to 
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The variance of the intensity fluctuations is obtained by setting ri = r2 = in eq. ( IB. 121 ). and making the change of variables, 
A6' = 0c - Oy, <d -{ex + dy)/2, X - r'j/rp and y - r'Jr^, to obtain 



{Af} + l^—\ dxdyxy I dA0 explixycos AO - 2A4x^y^(l + 2cos^ AO)] . 
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Thus the amplitude of the intensity variations is controlled only by the dimensionless parameter Arp. 
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The power spectrum of intensity fluctuations across the observer's plane is, 

W/(q) = J exp [-ir[ ■ q - 2A4rfq\\ + cos^C^,; - 0,)]} , (B.14) 

which is once again determined only by the dimensionless parameter Arp. 

No attempt is made to further analyse the characteristics of the intensity fluctuations in the regime of strong phase perturbations 
as this regime is unhkely to be encountered in the lensing of gravitational radiation even under the most extreme circumstances. 



Appendix C: Calculation of the temporal intensity fluctuation spectrum for lensing by stars 

The contribution to the temporal intensity fluctuation spectrum for oj/t; ^ Tp' from clustered stars is, 
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where we have written the stellar density profile to second-order in q, 1/^(^)1 ~ M (I - q R 15). The power spectrum is simplified 
by recognising that the last argument of the two hypergeometric functions is large for w < vr^^ . Using the asymptotic expansion, 
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keeping the leading order term in tolv, and neglecting terms of order « 1 we obtain, 

WeohM = -4rr(-r)^'sin(^)(r2M2/^V^r^l(^)', ^ < r;', 1/2 < y < 2.. 

For higher angular frequencies oj > vr^^, the temporal power spectrum is 
23-r 2^2^2 2^rr/_z\ 
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where, upon expanding the incomplete beta function for large z, 

Pz(«. 0) ~ — — 7T z (-Z) + — — — -zr(-z) . 

r(l -b) a + b-i 

neglecting terms OiR-^/r^) and keeping terms to lowest order in cj/v only one has 
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